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Abstract —The downlink of a wireless network where mnltl- 
antenna base stations (BSs) communicate with single-antenna 
mobile stations (MSs) using maximal ratio transmission (MRT) is 
considered here. The locations of BSs are modeled by a homoge¬ 
neous Poisson point process (PPP) and the channel gains between 
the multiple antennas of each BS and the single antenna of 
each MS are modeled as spatially arbitrarily correlated Rayleigh 
random variables. We first present novel closed-form expressions 
for the distribution of the power of the Interference resulting 
from the coexistence of one intended and one unintended MRT 
over the considered correlated fading channels. The derived 
expressions are then used to obtain closed-form expressions for 
the success probability and area spectral efficiency of the wireless 
communication network under investigation. Simulation results 
corroborate the validity of the presented expressions. A key result 
of this work is that the effect of spatial correlation on the network 
throughput may be contrasting depending on the density of BSs, 
the signal-to-interference-plus-noise ratio (SINR) level, and the 
background noise power. 

1. Introduction 

Multiple-input multiple-output (MIMO) systems have been 
instrumental in increasing the spectral efficiency and impro¬ 
ving the reliability of wireless links by combating channel 
fading and potentially reducing interference. Several multi¬ 
antenna techniques have been proposed, implemented, and 
standardized over the past decade, including diversity re¬ 
ception, beamforming, and spatial multiplexing. Diversity- 
oriented schemes enable multi-antenna transceivers to enhance 
the system performance by exploiting the diversity provided 
by the multiple channel fading variations. A classical antenna 
diversity technique is maximal ratio combining (MRC) HI, 
where the signals from the multiple receive antennas are 
weighted such that the signal-to-noise ratio (SNR) of their 
sum is maximized in the absence of interference, or when 
interference is treated as background noise. Maximum ratio 
transmission (MRT) m is the dual of MRC at the transmitter 
side, i.e. the transmit antenna weights are matched to channel 
fading, and exhibits the same performance. 

It has been long recognized that the theoretical performance 
gains of multi-antenna communication are degraded in practice 
due to correlated fading Ei-i]. Spatially correlated fading 
channels are usually encountered in multi-antenna systems 
employing not sufficiently wide separated antennas or with 
insufficient scattering around the transmitter. While spatial 


correlation has been considered as a drawback in both single- 
and multi-user MIMO systems (see e.g. |j9l), correlated fading 
may have a positive impact on spatial diversity techniques if 
opportunistic multi-user scheduling is employed Ho). 

The performance of diversity-oriented multi-antenna tech¬ 
niques in multi-cell networks has not been extensively studied 
until very recently, mainly due to significant difficulties in 
characterizing and modeling the co-channeFout-of-cell inter¬ 
ference. Recent advances in spatial network modeling using 
tools from stochastic geometry were key enablers for analyzing 
large spatial networks, departing from the classical simplifying 
assumption of Gaussian interference. Several recent papers 
have studied multi-antenna processing techniques in wireless 
networks with Poisson distributed interferers ini-llll, ig¬ 
noring, however, the effect of spatially correlated fading. In 
El, El, the effect of spatial interference correlation across 
diversity branches on the performance of single-input multiple- 
output (SIMO) links in Poisson networks was quantified. 
These works showed that multi-antenna diversity techniques 
suffer a diversity loss when spatial interference correlation is 
properly accounted for. Despite this progress, the performance 
characterization of multi-antenna techniques in spatially cor¬ 
related fading channels remains open, and is the main focus of 
this paper. In contrast to our paper, all aforementioned papers 
consider spatially uncorrelated Rayleigh fading, i.e. channel 
vectors are spatially white random vectors, and interference 
correlation occurs due to the common interference locations. 

In this work, we investigate the effect of spatially arbitrarily 
correlated Rayleigh fading on the performance of MRT in 
wireless Poisson networks, where interferers form a Poisson 
point process (PPP) and fading correlation at the antenna 
branches is present only at the transmitter. The latter con¬ 
sideration describes, for instance, a network set up in which 
multiple antennas are placed at a high-point base station (BS), 
while the mobile station (MS) is located in a rich scattering 
surrounding. We provide analytical expressions for the success 
probability and area spectral efficiency as a means to capture 
the effect of spatial fading correlation on the network perfor¬ 
mance. For that, we present novel closed-form expressions for 
the distribution of the squared norm of a normalized version 
of the dot product of two uncorrelated random vectors, each 
having correlated elements described by a common covariance 
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Fig. 1. The considered downlink wireless network. BS 0 communicates 
with the typical MS 0 under the presence of interference links (only N of the 
interference links are illustrated in this figure for exposition convenience). 


matrix. The validity of our theoretical analysis is corroborated 
through comparisons with equivalent computer simulations. 

Notation: Vectors and matrices are denoted by boldface 
lowercase letters and boldface capital letters, respectively. The 
transpose conjugate is denoted by j-)^, while I„ (n > 2) is the 
nxn identity matrix. | |a| | stands for the Euclidean norm of a, 
[a]ji represents a’s nth element, represents A’s (i,j)th 
element, and diagja} denotes a square diagonal matrix with 
a’s elements in its main diagonal. Nf, TZ and C represent the 
natural, real, and complex number sets, respectively, whereas 
E{-} is the expectation operator and P[-] represents probability. 
The Laplace transform (LT) of RV x is denoted as Cx{-} and 
X ~ CA/’(0,cr^) indicates that a; is a circularly-symmetric 
complex Gaussian RV with zero mean and variance 
Finally, r(-) is the Gamma function ifTTl eq. (8.310/1)], Ei(-) 
is the exponential integral Elleq. (8.211/1)1, and 2^1 
denotes the Gauss Hypergeometric function iflTl eq. (9.14)]. 

II. System and Channel Model 

We consider the downlink of a wireless network, as shown 
in Fig. [T] where BSs are arranged according to a spatial 
homogeneous PPP $ C TZ? of density A in the Euclidean 
plane. Each BS is equipped with nx > 2 transmit antennas and 
communicates with one single-antenna MS at a fixed distance. 
Our results can be easily extended for random link distances, 
e.g. assuming a geographically nearest BS connectivity model. 
All simultaneous downlink transmissions are assumed per¬ 
fectly synchronized, and each BS k, with k G Af+, processes 
individually its complex symbols Sk with a linear precoding 
vector ffc G C"'^ before transmitting. It is also assumed that 
for each value of k, ||ffc|| = 1 and E{|sfcp} < P, where P 
denotes the transmit power of each BS. 

Let G be the narrowband channel vector for the 
wireless link between MS k and BS £ and also let denote 
the Euclidean distance between these nodes. The baseband 
received signal at MS k can be mathematically expressed as 

Vk = (1) 

where the standard power-law pathloss model l{r) = with 
common pathloss exponent a > 2 for all links is assumed. 


Ik is the cumulative interference from all BSs except the 
associated BS k, which is given by 

Ik= ^ (2) 


and rife ~ CAf{0, a^) denotes the additive white Gaussian 
noise (AWGN). Based on the considered system model de¬ 
scribed by O and the signal-to-interference-plus-noise 
ratio (SINR) at MS k can be computed as 


-1 


7fc = »'fe,fe|hfc,fcffe| 




c \2 


TT (3) 


The elements of each channel gain vector are assumed to be 
arbitrarily correlated and in particular, each hfe^^ is modeled as 
= R^/^hfe_^, where hfe^^ G is a standard circularly- 
symmetric jointly-Gaussian random vector Ha Sec. 3.7] and 
R = E{hfe_fh^^} denotes the common covariance matrix of 
hfe^^’s. In addition, we assume that R is real symmetric having 
ones in its main diagonal as well as that hfe/’s are independent 
across the different BSs and of both rj. ^’s and rife’s. 

In this paper, we assume that each BS k knows perfectly 
hfe fe and treats interference as background noise. BSs do 
not cooperate 01 and each of them acts selfishly aiming at 
maximizing its own SNR. In order to accomplish that, each 
BS k is assumed to perform channel-matched precoding, a.k.a. 
MRT. Hence, ffe at each BS k is designed as |l2] Sec. Ill] 



HI. Pereormance Analysis 


( 4 ) 


In this section we present analytical expressions for the suc¬ 
cess probability and the area spectral efficiency of the conside¬ 
red downlink wireless network. The network performance is 
analyzed assuming a typical MS 0 located at the origin. As 
a result of Palm probabilities (see li20ll and references therein 
for details) and the stationarity of the PPP, the statistics of 
signal reception at the typical MS 0 is seen by any MS in the 
considered wireless network. 


A. Success Probability 

The success probability of a typical MS is defined as 

Rjuc = p [70 > 7th] ( 5 ) 


where 7 th denotes a predetermined SINR threshold. By sub¬ 
stituting (|4|i into (O for MS 0 and then into (|5]l, Pguc can be 
expressed using the cumulative distribution function (CDF) of 
II hop IP and the expectation of the real RV that represents 
the aggregate interference, and is given by 

'^'0= E (6) 

ie$/{o} 


with the complex RV gi defined as 

Th u. 

IlhMil ■ 


( 7 ) 















In particular, Pguc can be obtained as 


-fsuc - 1 


^Ilho.oP 


^0.07th 



( 8 ) 


where E^j,{-} represents the expectation of 'I'o- By setting the 
Nakagami parameter m = 1 in [3 eq. (20)] (resulting in the 
Rayleigh case) and integrating, a closed-form expression for 
the CDF of 11 hop IP is easily derived as 


^iiho,oip(®) - x!n I ^ 

t=l i=l 


Ml 

[^Al 


-1 


1 — exp 



(9) 

where fi G 72.”'^ contains the distinct eigenvalues of R in 
decreasing order (the non-diagonal elements of R lie in [0,1)). 
Hence, plugging @ into ([8]l and using the dehnition of the 
LT of ^'o, the following expression for Psuc can be deduced 
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To obtain a closed-form expression for Psuc based on the 
latter expression (fTOb , we hrst derive an exact as well as an 
approximate analytical expression for the probability density 
function (PDF) of V i € J\f+ by means of the following 
two theorems. 


Theorem 1. The PDF of the RV g = jg^p V i G A/+ for the 
special case where np = 2 is given by 

fgix) = ^-^[Pil+P,x) - P{1-p,x)] (11) 

where p = [R]j^ ^ C [0,1) and function P{p,x), with p > 0 
and X > 0, is dehned as 


P(?7,x) 


V 

4(2 - p) 


exp 






Proof: The proof is provided in Appendix 


( 12 ) 


Theorem 2. The PDF of g for the general case np > 2 can 
be approximated by the following exponential PDF 


/g(x) = — exp- 


where the scale parameter cTg is given by 


= E 


Mi 

^1 Mn ■ 


(13) 


(14) 


Proof: See Appendix iBl ■ 

Remark 1: Theorem |2] also holds for the special case where 
R = I„ (spatially uncorrelated Rayleigh fading), i.e. [/r]^ = 1 
V n = 2, 3,..., nx (identical eigenvalues). For this special 



Fig. 2. Histogram of the positive real RV g for the special case of nx = 2 
and vaiious values of the correlation coefficient p. 



Fig. 3. Histogram of the positive real RV g for a correlation coefficient 
p = 0.8 and vaiious values of nx- 


case, (fT4li simplihes to Cg = 1 and Theorem |2] provides 
the exact PDF of g for nx > 2, which coincides with the 
equivalent result in ifT^ and ||2T1 . 

The validity of the exact PDF expression in Theorem [T] and 
the tightness of the approximate PDF expression in Theorem|2] 
are demonstrated in Figs.|2]and[3 respectively. In these hgures 
the histogram of g is plotted for various exponential correlation 
matrices, i.e. [R], = p\^-3\ a, and various values for nx- As 

clearly shown in Fig. |2] the curves for the exact PDF expres¬ 
sion for nx = 2 given by (fTTTl match perfectly with the curves 
for the equivalent empirical PDF. In addition, as depicted in 
Fig. [3l the curves for the approximate PDF expression (foT l 
for nx > 2 exhibit a good ht with the empirical ones. This 
ht improves as p and/or nx increase (or equivalently as the 
histogram becomes more heavy tailed). The latter substantiates 
that Theorem |2] provides an accurate approximation for the 
PDF of g for multi-antenna communication networks with 
increased spatial fading correlation and very large numbers 
of antennas (massive MIMO). 

Since Ipip’s are independent and identically distributed RVs 






























and independent from the PPP $, the LT that appears in ( fTOl i 
can be obtained using similar steps to ll^ Theorem 1] as 


(15) 


pCO 

(s) = exp —27rA / TL^s, x)fg{x)dx 
L Jo 

where function 'H{s,x) is defined as 

poo 

'H{s,x)= [l-exp{-sxy~°‘)]ydy 
Jo 


(i) (sx) 
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-n 1-- 
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(16) 


In (fThl l. (i) follows from the change of variables y ^ u and 
after some algebraic and calculus manipulations. Substituting 
(O into ( fT5l l. then replacing either (fTTT i or (fT3l l. and finally 
using the definition of the fractional moment, a closed-form 
expression for the LT of 'I'o can be obtained as 

(s) = exp (17) 

where the parameter ^ > 0 is given by 

e = ttT (^1 - E{g 2 /“}. (18) 


A closed-form exact expression for the special case of tt-t = 2 
and a closed-form analytical approximation for the general 
case riT > 2 for the fractional moment E{g^/“} that appears 
in (fTSl l are presented in the following two corollaries. 


Corollary 1. A closed-form exact expression for E{g^/“} for 
the special case where nx = 2 is given by 

IE{g 2 /a} ^ ^ P _ g (X - p, a)] (19) 

2p 

where function Q{r],a) is defined as 


Q{v,a) = 


2/a-l-l 
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( 20 ) 


Proof: Starting from the definition of the expectation 
operator and using (fTTIi and (fT2li . (l20l i can be easily obtained 
with the use of ifTTl eq. (3.381/4)], ifTTI eq. (6.228/2)] and after 
some elementary algebraic manipulations. ■ 


Corollary 2. A closed-form approximate expression for 
E{g^/“} for the case where nx > 2 is obtained as 

E{g 2 /“} = a 2 /“r + 1 ) ■ ( 21 ) 

Proof: Substituting (fT3l l and (fT4l i into the definition of 
the expectation operator and making use of ini eq. (3.381/4)] 
completes the proof. ■ 

Remark 2: Similar to the remark for Theorem |2] for the 
special case where R = I„ V n = 2, 3,..., rtx, Corollary |2] 
presents an exact expression for E{g^/“} by setting Ug = 1 . 


To derive a closed-form exact expression for Psuc for the 
special case nx = 2 , (fTTl i is substituted into (fTOl) yielding 


P.„=^i:iLL^exp. 




2p 


X exp < —^Arp 


0,0 


^ 0^,0 7th 

P[l-(-l)V] 

-| 2/q; ' 
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( 22 ) 


where ^ is given by (HI after substituting (O with 
Similarly, substituting (fTTl) using (fTsT i and (l2Tll into (fTol l. the 
following closed-form approximation for Pguc is obtained 

riT ttt / r 1 \ —1 

[P']^ 




l=\ i=l 


1 — exp — 


7th 
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exp 




2/0 


(23) 


B. Area Spectral Efficiency 

Using the previously derived expressions for Psuc, the 
area spectral efficiency of the considered downlink wireless 
network is straightforwardly obtained as 


ASE = Alog2(l + 7th)Psuc. (24) 

By substituting (l22l) for nx = 2 and ( 12^ for nx > 2 in (l24l) . a 
closed-form exact and a closed-form approximate expression, 
respectively, for ASE can be easily derived. 


IV. Performance Evaluation Results 

In this section, numerically evaluated results for the perfor¬ 
mance expressions provided in Section [111] are presented and 
compared with equivalent results obtained by means of com¬ 
puter simulations. In particular, we numerically evaluate both 
the exact and the approximate expressions for Psuc given by 
(I 22 I 1 and ([23), respectively, as well as those for ASE obtained 
from (i24l) . Without loss of generality, in the performance 
evaluation results that follow, exponential correlation matrices 
have been considered and it has been assumed that a = 3.5, 
ro_o = 1 and cr^/P = 1 . 

In Fig. H is plotted versus 7 th for A = 10 nx = 2 
and 5, and for p — {0.01, 0.7, 0.95}. As shown for all cases, 
the numerically evaluated results for the derived exact (nx = 
2) and the approximate (nx = 5) Psuc match perfectly and are 
sufficiently close, respectively, with equivalent simulations. In 
addition, it is demonstrated that Pguc improves with decreasing 
7 th and increasing nx. Interestingly, there exists a critical 7 th 
value, 7 *j^, (crossing point in the figure) that increases with 
increasing nx that determines the impact of p on Psuc- More 
specifically, for 7 th < 7 th, Psuc improves with decreasing p. 
However, when 7 th > 7 ^ 1 , p decreases Pguc degrades. 

By setting 7 th = 3 dB in Fig. [Sj ASE is illustrated as 
a function of A for nx = 2 and p = {0.01,0.7,0.95}. As 
shown in this figure, there exists a perfect match between the 
evaluated results for the derived exact ASE and equivalent 
simulations. In addition, the derived approximation appears to 





















Fig. 4. Success probability, Psuc, as a function of the SINK threshold 7 tii 
in dB for A = 10“^, nx = 2 and 5, as well as for p = {0.01, 0.7, 0.95}. 



Fig. 6. Area spectral efficiency, ASE, as a function of the BS density A for 
7 tli = 5 dB, nx = 5 and 20, as well as for p = 0.7 and 0.95. 



resulting from the coexistence of one intended and one un¬ 
intended MRT over the considered fading channels were first 
presented. Based on these expressions and using tools from 
stochastic geometry, we further derived closed-form expres¬ 
sions for the success probability and area spectral efficiency 
of the downlink network under investigation. Our theoretical 
analysis, which is corroborated with simulation results, shed 
light on how spatial fading correlation affects the network 
performance depending on the interplay between the density 
of BSs and the target SINR threshold. Interestingly, depending 
on the background noise level, spatial fading correlation may 
increase or decrease the network throughput, having a phase 
transition for a certain SINR threshold value. 


Fig. 5. Area spectral efficiency, ASE, as a function of the BS density A for 
7th = 3 dB, nx = 2, and p = {0.01, 0.7, 0.95}. 


be a tight lower bound, which tightness increases as p tends 
to either of its two extreme values (i.e. p —>• 0 or p —i 1) as 
well as when A decreases. It is also demonstrated that ASE 
degrades with increasing p and its maximum value happens for 
lower values of A as p increases. The same performance metric 
over A is illustrated in Fig.|6]for cases where p = 0.7 and 0.95, 
n-T takes the values 5 and 20, and 7 th is increased to 5 dB. 
In this figure, it is shown that ASE improves with increasing 
riT- In contrast to Fig. |5l as p increases, ASE is shown to 
improve. Moreover, the maximum value for ASE happens for 
lower values of A as p decreases. As also discussed in the 
description of Fig. |4l the latter behavior implies that, for the 
parameters’ setting in Figs. |5] and |6] p has different impact on 
ASE. Finally, it is noted that, as demonstrated in Fig. |6l the 
impact of p on ASE increases with increasing ut- 

V. Conclusion 

In this paper, we have analyzed the performance of MRT 
in a wireless Poisson network assuming spatially arbitrarily 
correlated Rayleigh fading channels. Novel closed-form ex¬ 
pressions for the distribution of the power of the interference 


Appendix A 
Proof of Theorem[T] 

Starting from © and using the definition of the channel 
gain vectors ho^ and i in Section HH pi V { G A+ can be 
equivalently expressed as gi = with the random vector 

Wi G defined as 


Wi = 



R 


1/2 




(A.l) 


It is obvious that = w^ho,i V * G A/+ conditioned on hpi is 
a linear functional of the standard circularly-symmetric jointly- 
Gaussian random vector ho,i, and hence it is a circularly- 
symmetric Gaussian RV m Sec. 3.7]. The mean value of 
each Pi can be easily verified to be independent of ^ and 
given by 

nx 

V- = ^E{[ho,.]J[wf]J:^^0 (A.2) 

n—1 


where (b) follows from our assumption on the expectation of 
the elements of in Section |II] In addition, it can be shown 
that each |pi|^ conditioned on is exponentially distributed 
with scale parameter obtained as 


^{\9^f \K^} 


RhM ■ 


(A.3) 


































In the derivation of the latter expression, we have used the 
notation = RR and the fact that E{ho,ih^j} = 12- 
To obtain an analytical expression for the distribution of 
\gi\^ V i G A/+, we need first to derive an analytical expression 
for the distribution of the generalized Rayleigh quotient that 
appears in the right-hand side of (IA.3T i.e. for the positive 
real-valued RV n = hpjR^hi_i/hpjRhi,i. Using the results 
of II 22 I . a closed-form expression for the PDF of each r.^ can 
be obtained as 


fri (y) 


(2-J/)^ l-p<y<l + p 
0, otherwise 


It is noted that, as shown in (IA.4b . the PDF of each depends 
only on p, and hence is identical for all r^’s (by assumption, 
R is common for ^ V i G A/+). By using the exponential 
conditional PDF of \gi^ and (IA.4b . the unconditional PDF of 
each \gi^ can be finally computed as 


manipulations yields (fTTT i. 


f 1 f 

x\ 

= - exp 

— 

Jo y V 

yj 

l_p 2 |.i+p exp 

Ji-p 

2/(2 

eqs. (3.352)] 

and 


(a) 


(A.5) 


dy, 


Appendix B 
Proof of Theorem[2] 

Both R^ and R that appear in the scale parameter described 
by (IA.3b are real symmetric matrices, and hence their eigen¬ 
value decomposition is computed as R^ = Vdiag {v>} and 
R = Vdiag{/x}V^, respectively, with V G being 

an orthogonal matrix and [i/]^ = [/x]^ V n = 2, 3,..., nx- By 
introducing the random vector i, which can be ea¬ 

sily shown to be circularly-symmetric jointly-Gaussian, (lA3l l 
can be re-expressed as I = E{|pi|^ \ ri\ = ri, 

where the positive real-valued RV that was introduced in 
the previous appendix is hereinafter equivalently given by 




M.|[q4 

Mn l[q*]r 


(B.l) 


To derive an approximate expression for the PDF of each 
|pi|^, we make use of one of the results presented in Ap¬ 
pendix according to which, the PDF of each is an 
exponential one with scale parameter depending on ha, or 
equivalently on ri as expressed as (iBdT i. We propose the 
approximation of the PDF of V x G N+ with an exponen¬ 
tial PDF having as a scale parameter the E{|pi|^} = E{ri}, 
which does not depend on ri. To this end, E{|( 7 i|^} can be 
approximated as follows 


E{|p,| 2}^E £[/x] 

ln=l 


= £[/x] 


^ n—1 



(B.2) 


where uji^n = |[qi]n,|^- By using the fact that E{a;i_„} = 1 V 
i G A/+ and V n = 2, 3,..., tit, the right-hand side of (IB. 21) 
can be expressed V 1 G as in (fT4l i. 
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